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ABSTRACT 

We investigate whether a correlation exists between the concentration (200) and spin 
(A) parameters of dark matter haloes forming in ACDM N-body simulations. In partic- 
ular, we focus on haloes with virial masses in the range 10 11 ^ M v i r //i _1 M Q ^ 2x 10 12 , 
characteristic of the hosts of Low Surface Brightness (LSB) and High Surface Bright- 
ness (HSB) galaxies. Our analysis demonstrates that 200 and A are anti-correlated. If 
we assume that a galaxy disk forms in one of these haloes from baryons that approx- 
imately conserved angular momentum during their dissipation into the rotationally 
supported disk, then it is possible to estimate the disk's central surface density. For 
an appropriate choice of mass-to-light ratio, we can place constraints on the central 
surface brightness of the disk and hence identify the analogues of the host haloes of 
LSB and HSB galaxies. We find that our LSB galaxy analogues occupy haloes that have 
lower concentrations than might be expected based on consideration of the M v i r -c v i r 
relation for the ACDM cosmology. The distribution of concentrations peaks at 200~ 6, 
in good agreement with observational data, although there are important differences 
between the shapes of the simulated and observationally inferred distributions. This 
suggests that LSB galaxies inhabit a biased subsample of the halo population, and that 
this bias may be an important ingredient for resolving the current debate regarding 
the structure of LSB galaxy dark matter haloes. 

Key words: galaxies: haloes - cosmology:theory, dark matter, gravitation - methods: 
numerical, N-body simulation 



1 INTRODUCTION 

The Cold Dark Matter (CDM) model is widely accepted 
as the paradigm within which to understand the forma- 
tion of structure in the Universe. A defining characteris- 
tic of this model on the scale of galaxies is the predic- 
tion that central densities of dark matter haloes rise as 
p oc r a , with a « — 1, over the radii which govern the in- 
ner rotation curves of galaxies. Over the last decade, nu- 
merous studies - principally numerical - have sought to 
quant ify and understand this d i vergen t or cusyy behaviour 
(e.g. INavarro. Frenk fc White! 1 19961: iMoore et aD Il998t 
iNavarro et al.ll2004t) . In particular, the continued assertion 
that the rotation curves of dark matter dominated sys- 
tems such as low surface brightness (LSB) galaxies require 
mass profiles with finite constant density cores or shallow 
"cusps" (e.g.lFlores fc Primacklll99l lMoorelll994l: iBurkerd 
1993:lde Blok fe McGaugrJl997l:TMcGaugh fc de Blokll998t 
de Blok et al]l2003h has prompted many detailed investi- 
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gations into the consequences of these observations for the 
CDM model. 

This disparity between the predicted and inferred 
mass profiles appears to call into question the most basic 
assumption of the CDM model, namely the nature of 
the dark matter, and constitutes a critical challenge to 
the paradigm. Indeed, reconciling theory and observation 
has come to represent a central problem in contemporary 
cosmology and is arguably one of its most contentious 
issues, as illustrated by recent analyses of the systematic 
effects that could affect high resolution rotation curve 
data of LSB g alaxies and the fittin g procedures used in 
mass modeling. Ide Blok et alJ (l2003t hereafter DB03) have 
argued that the signature of CDM cusps should be clear 
in the data and so the model's predicti o ns are inconsistent 
with observations, while lHavashi et aD J2004allbh . in some 
cases using overlapping data, find no compelling evidence 
for disagreement between cuspy CDM haloes and observed 
rotation curves. These studies focused on the shape of LSB 
galaxy rotation curves, which rise approximately linearly 
with radius in the manner one would expect if the mass 
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density profile contained only a shallow central cusp or a 
constant density core; if the underlying profile is strongly 
cuspy, it is argued, the rotation curve should rise more 
rapidly at small radii. However, the rate at which the 
rotation velocity increases with radius is sensitive to how 
centrally concentrated the mass distribution within the halo 
is - the more centrally concentrated the halo, the more 
rapidly its rotation curve will rise. Thus, a slowly rising 
rotation curve could be interpreted as indicative of a halo 
that is weakly centrally concentrated. 

The seminal work of iNavarro. Frenk fc White! l)l99fit 
hereafter NFW) formalised this concept of concentration. 
Based on a comprehensive suite of simulations, they estab- 
lished their universal mass profile that could describe the 
structure of an average dark matter halo in dynamical equi- 
librium, independent of its mass or the choice of cosmological 
parameters and initial density perturbation spectrum. The 
circular velocity of a NFW halo can be expressed as 

Vc(x) _ j 1 ln(l + C vir x) - C vil x/(1 + CyirZ) 
Kir ~ V X M 1 + C ™) ~ <W(1 + C vir ) ' 

where x = r/r v i r is the radius r normalised by the virial 
radius r v - lT , c v i r = r v i r /r s is the concentration parameter 
where r s is the scale radius, and Vnt = \J GM v i r /r v i T is 
the virial velocity of the halo. From equation Q we observe 
that the circular velocity rises as r 1 ' 2 at small radii, and the 
larger c v i r , the greater V c (x)/V v i r - 

NFW found that, on average, concentration decreased 
with increasing virial mass, and they understood that this 
reflected the mean density of the Universe at the time 
of the halo's formation; more massive systems collapse at 
later times when the mean den sity of the Universe is lower. 
Navarro. Frenk fc White Jl997f) and subs equent studies by 
Bullock et al.1 fcOOlal. hereafter B01) and lEke et ail feOOll 
hereafter ENS) formulated prescriptions that allow the con- 
centration of a given halo of a given mass in a given cosmol- 
ogy to be predicted. Within the B01 or ENS formulations, 
we expect a typical galaxy mass halo (2 x 1O 12 M ) forming 
in the ACDM cosmology to have a concentration parameter 
c vir ~ 10. 

If LSB galaxies occupy dark matter haloes of the kind 
that form in the ACDM model, the n we might expect typi - 
cal concentrations of c v j r ~ 10. Both lMcGaugh et alJ <l2003t) 
and DB03 have considered the distribution of concentrations 
obtained from fits to their LSB galaxy rotation c urve data, 
and as illustrated in Figure 3 of lMcGaueh et all (I2003T) and 
Figure 1 of DB03, the distribution shows two well defined 
peaks, one at c v j r ~ 5, the other at c v i r = 0. The preference 
for the data to favour such a low value of concentration ap- 
pears troubling for the CDM model. The mass-concentration 
relation, predicts that c v ir ~ 5 is typical of low mass clusters 
with Mvir ~ 1O 14 M0 rather than low mass galaxies, and the 
1 — a scatter in concentration at a given mass is expected 
to be 50% at most. 

Implicit in this comparison is the untested assumption 

1 The virial radius defines the extent of a dark matter halo and is 
such that it encloses an overdensity A v j r times the critical density 
of the Universe at a given epoch; for the ACDM model, A v i r ~ 
100. 



that the haloes of LSB galaxies are typical for their mass. 
LSB galaxies are considered ideal for probing the structure 
of dark matter haloes because the baryons contribute mini- 
mally to the rotation curve, greatly facilitating the process 
of decomposing the rotation curve into baryonic and halo 
contributions. However, if the haloes of LSB galaxies are bi- 
ased with respect to the general halo population, then this is 
a biased comparison. Indeed, it has been noted observation- 
ally that detailed mass models of LSB galaxies reveal dark 
matter haloes that are less dense and more extended than 
those of high surface bright ness (HSB) counterp arts of simi- 
lar masses (see the review of lBothun et alll997f) . suggesting 
lower concentrations than might be expected. Furthermore, 
there are theoretical grounds to expect that LSB galaxies in- 
habit a special class of haloes. The low surface brightness of 
LSB galaxies is a consequence of their low central density. 
If we make the reasonable assumption that specific angu- 
lar momentum of the baryons is approximately conserved 
during their dissipation into a rotationally supported disk, 
then the scale length of the resulting disk is related to the 
angular momentum of the dark matter halo, and we antici- 
pate that LSB galaxies will occupy dark matter haloes with 
larger spin parameters (a measure of the degree to which the 
halo is supported by angular momentum; see §[2J) than their 
HSB counterparts. These considerations suggest that LSB 
galaxies occupy dark matter haloes with low concentrations 
and high spin p arameters, echoing the conclusions of previ- 
ous studies (e.g. IDalcan ton et alJll997t Ijiminez et al. 1998; 
iMo. Mao fc Whitdll99Sft ~ 

In this paper, we investigate whether the correla- 
tion between concen t ration and spin, as e nvisag ed by 
e.g. iDalcanto n et al 1 lll997tl . | Jiminez et al] Jl998), and 
IMo. Mao fc White! ijlflflaT is sufficie nt to explain the distri- 
bution of concentrations observed bv lMcOaugh et ail ( 2003) 
and DB03 for their sample of LSB galaxies. For this purpose, 
we have analysed a statistical sample of haloes, with masses 
in the range 10 11 ^ M„ ir //i _1 M < 2 x 10 12 , forming in 
a high resolution cosmological ACDM N-body simulation, 
and observe a trend for a halo's spin parameter to increase 
with decreasing concentration. Using the a nalytic prescrip- 
tion presented in lMo. Mao fc White! dl998l) . we use the halo 
spins to estimate disk scale lengths and surface densities for 
our hypothetical sample of galactic disks, and hence con- 
struct the distribution of concentrations we might expect 
for the host haloes of LSB galaxies. Although there are im- 
portant difference s between the distribu tion we recover and 
those presented in lMcGaueh et ail fcOOSt) and DB03, we ob- 
serve that the peak at c v i r ~ 5 is reproduced in our distri- 
bution. We present the results of our analysis in the next 
section (§0 and discuss their significance in §0 



2 RESULTS 

Our results are based on analysis of the z~0 output of a 
512 3 particle cosmological N-body simulation of a periodic 
volume, 50 h -1 Mpc on a side, that was run assuming the 
ACDM cosmology (Q = 0.3, fl A = 0.7, tt b = 0.045, h = 0.7, 
a s = 0.9) using the GADGET2 code (Springel et al., in 
preparation). A (comoving) gravitational softening length 
of 5/i _1 kpc (approximately l/20th the mean interparticle 
separation) was adopted, and the chosen parameters imply 
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Figure 1. The correlation between spin parameter, A, and con- 
centration, C200. We include all haloes with masses in the range 
10 11 < M vil /h- 1 M (3 ^ 2 X 10 12 . See text for further details. 



a particle mass of 7.75 x 10 7 h -1 Mq. The mass resolution 
of the simulation is sufficient to resolve a 2 x 10 12 fo -1 Mq 
halo with ~ 26, 000 particles. Haloes are i dentifi ed using a 
friends-of-friends (FOF; e.g. iDavis et alj |l985) algorithm 
with a linking length of b — 0.2 times the mean interparticle 
separation. We not e that this simulation form ed the basis 
of recent studies bv lBailin fc Steinmetzl J2004al fbh. 

In order to facilitate comparison to observations, we 
define the concentration with respect to the radius that en- 
closes a spherical overdensity of 200 times the critical den- 
sity (as in the original NFW papers), rather than the cos- 
mology dependent A ~ 100. This has the effect of lower- 
ing the measured concentration of the halo by ~ 30% (ap- 
proximately independent of concentration) relative to that 
predicted by either the B01 or ENS prescriptions, reflect- 
ing the smaller measured virial radius relative to a fixed 
scale radius. Concentrations are determined by identifying 
best fitting parameters (r*2oo,C2oo) for the enclosed mass at 
(0.06,0.12,0.25,0.4,0.6,1.0) r vir ; we find that this method is 
efficient and it produces accurate and robust estimates of 

C200- 

The standard definition of the spin parameter, 



A 



J|£| 1/2 



(2) 



GM 5 / 2 ' 

where M, J and E are the mass, total angular momentum 
and energy of the halo respectively, depends on binding en- 
ergy, which is both a slow operation (~ iV 2 ) and ill defined 
in the context of a cos mological simul a tion. In stead, we have 
chosen to compute the lBullock et al . (2001b) measure, 

A' = -J- , (3) 

V2MVR 

where V is the circular velocity at radius R; this can be 
corrected by a factor A = A'/(c2oo) 1//2 , 



/(c2oo) = 2/3 + (c2oo/21.5) a7 

iMo. Mao fc Whitelll99dlBullock et al]l2001bl) . 



(4) 



centrations C200 relative to their spin parameters A for all 
FOF groups with masses between 10 11 and 2 x 10 12 h -1 M , 
as shown in Figure Q Although there is large variation in A 
at a given C200 , we note a discernible trend for A to increase 
with decreasing C200; the median increases from A mc d — 0.03 
at C200 = 10 to A mc d — 0.05 at C200 = 5. This is pre- 
cisely the behaviour we would expect the host haloes of 
LSB gala xies to display, in agreement with pr evious work 
(e.g. iDalcanton et alill997t iJiminez et aljll998T) . 

Having established a correlation between A and C200, we 
wish to derive the quantity that characterises a LSB galaxy, 
namely its central s urface brightness mo j our p articular treat- 
ment follows that of lMo. Mao fc White! (I1998T) . If we assume 
that the angular momentum of the baryons is approximately 
conserved during their dissipation into a rotationally sup- 
ported disk, we can estimate the angular momentum of the 
disk Jd and hence its extent for a given mass Md- We as- 
sume that the disk's mass and angular momentum are fixed 
fractions of those of the host halo's mass and angular mo- 
mentum, and choose j d = mj = 0.05 such that Md = raid M 
and Jd = jd J- Accordingly, the central surface density £0 
is given by 



E 



where Td is the scale length of the disk, 

»"200 / jd \ , f, s-l/2 £ 1 \ ■ \ 

rd = —^\ A/ C200 fR(X,C2oo,m d ,jd), 

V2 \m d J 



where 



/fl(A,C2oo,?rid,id) = 2 



2 Vc(r d u) 
V200 



du 



(5) 



(6) 



(7) 



and we use the fitting formula of lMo. Mao fc White! (Il99gft 
(see their equation (32)). The central surface brightness mo 
can be estimated by assuming a particular choice of local 
mass-to- light (M/L) to t ratio, which is a combination of the 
gas {M/L)m and star (M /L)* ratios. Ide Blok et al] (Il997l) 
indicates a local (M/L)hi ~ 1, and de Blok (2004, private 
communication) suggests (M/L), ~ 1; in what follows, we 
adopt (M/L)tot = 2. 

By definition, a LSB galaxy is one whose central 
surface brightness satisfies ^ 23 mag arcsec -2 ; typical 
B-band values lie in the range 23 /10 ^ 24 mag arcsec -2 
for LSB galaxies, compared with 21 ^ mo ^ 22 mag 
arcsec -2 for HSB galaxies. Converting to central surface 
densities, we select LSB and HSB galaxies by applying cuts 
of 1.5 ^ log Eo ^ 1.9 and 2.3 ^ logEo ^ 2.7 respectively, 
where Eo is expressed in units of Mq pc -2 . 

In Figure [5] we demonstrate how central surface den- 
sity Eo correlates with concentration C200, where Eo is eval- 
uated using equations ijjjj and ©. As expected, there is 
a clear trend for haloes with lower concentrations to host 
disks that have lower central surface densities; the median 
increases from Eo, mc d — 10 2 Mq pc -2 at C200 = 5 to 
E , mc d — 10 3 Mq pc -2 at C200 = 10. We quantify this trend 
by computing the mean value of C200 in bins of width 0.5 
in log Eo and performing a linear regression over the range 
1 ^ log Eo ^ 3.5 to obtain the relation, 



We begin by considering the distribution of halo con- 



C200 



3.4 + 2.1 log E , 



(8) 
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Figure 2. Correlation between central surface density, So, and 
concentration, C200- The solid line represents the variation of the 
mean value of C200 with So- As in Figure we consider only 
haloes with masses 10 11 Mvir/'i~ 1 M0 < 2 x 10 12 . See text for 
further details. 



0.3 



o 0.2 



0. 1 



0.0 













Blok e 


ol. (2003)- 












LSB halo 

HSB halo 


somple 
sample 





















5 10 15 20 



Figure 3. Distribution of concentrations. See text for further 
details. 

shown by the solid curve in Figure [5] outside of this range, 
C200 ~ 5 for log Eo < 1 and C200 ~ 10 for log Eo ^ 3.5. This 
correlation is due in equal parts to the direct dependence 
of Td on C200 and on the indirect dependence via the spin 
parameter A. For a fixed Eo, we find the distribution of C200 
has a width of ~ 2.5 — 3, essentially independent of the 
precise value of Eo. 

Finally, in Figure [3] we show distributions of concen- 
trations corresponding to the host haloes of LSB (red his- 
togram; hereafter the LSB distribution) and HSB (blue his- 
togram; hereafter the HSB distribution) galaxies by apply- 
ing the observationally motivated cuts in central surface 
density (1.5 < logE < 1.9, 2.3 < logE < 2.7 for LSB 
and HSB disks respectively) discussed above. For compari- 
son, we show also the distribution derived from fits to high 
resolution LSB galaxy rotation curves presented in Figure 1 
of DB03 (hereafter the DB03 distribution). 

There are two main points worthy of note in this figure. 



Firstly, we find that the LSB and HSB distributions peak 
at different values of C200; the peak value of the HSB distri- 
bution occurs at C200 ~ 8.5, compared with C200 ~ 6 for the 
LSB distribution. Secondly, the values of C200 at which the 
LSB and DB03 distributions peak are in remarkably good 
agreement - C200 ~ 6 - although we note that the overall 
shapes of the respective distributions differ; the DB03 distri- 
bution has an additional well defined peak at C200 = and a 
tail extending to C200 = 20. However, the peak at C200 = in 
the DB03 distribution is an artifact of the fitting procedure; 
when the NFW fits were clearly unphysical, such as when 
they preferred negativ e values of C200, the valu e of C200 was 
arbitrarily set to 0.1 ide Blok fc Bosmal I2OO2I . Therefore, 
these haloes should be characterized as haloes which cannot 
be fit by the NFW form, rather than haloes with a well- 
defined concentration of 0.1. We also note that the M v i r -C2oo 
relation implies that very high concentration haloes are ex- 
clusively very low mass systems. In particular, the observed 
halo with C200 ~ 19 is likely a dwarf galaxy whos e mass 
falls well below our low mass cutoff of 10 11 ft" 1 M Q iPowerl 
12003ft . 

We can estimate the likelihood that the DB03, LSB 
and HSB samples are drawn from the same distribution by 
applying the K-S test. If we include the entire DB03 sample, 
we find that it is inconsistent with both the LSB (Pks = 
0.04) and HSB (2 x 10" 6 ) samples; ignoring the peak at 
C200 = 0, we find reasonable agreement with the LSB sample 
(0.11) but not the HSB sample (0.002); and finally, ignoring 
both the peak at C200 = as well as the high-C2oo tail, we 
find improved agreement with the LSB sample (0.15) but 
similarly poor agreement with the HSB sample (0.001). 



3 DISCUSSION 

The status of LSB galaxies as dark matter dominated 
systems is a well established observational result (e.g. 
Bot hun et £01119971) and their rotation curves can provide, 
in principle, clean probes of the distribution of dark mat- 
ter within their host haloes. However, as the divergent 
conclusions reached by the recent studies of DB03 and 
lHavashi et alJ i2004al) demonstrate, the question of whether 
the shape of LSB galaxy rotation curves can be used to 
place strong constraints on the nature of the dark matter 
and hence the validity of the ACDM model on small scales 
continues to be a matter of heated debate. We have not at- 
tempted to address this particular issue in our study - the 
mass resolution of our simulation is insufficient to reliably 
recover halo density profiles on the scales that are of most in - 
terest for comparison with observations dPower et al.ll20(53) 
- but it is unlikely that their structure will differ significantly 
from that predicted by NFW and more recent studies (e.g. 
iNavarro et al.ll2004l) . Instead we have focused on the related 
topic of the degree to which haloes of the kind predicted 
to host LSB galaxies in the ACDM cosmology are centrally 
concentrated. 

Our analysis indicates that a correlation exists between 
concentration C200 and spin A for haloes with masses in the 
range we have studied, 10 11 < Mvi r //i _1 M© < 2 x 10 12 , 
such that haloes with low values of C200 tend to have high 
values of A, and vice versa. By assuming that the baryons 
associated with these I0W-C200 , high-A haloes approximately 
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conserve angular momentum during the formation of the 
galaxy disk, we recast this trend as one between the central 
surface density of the disk, Eo, and the concentration of its 
host h alo, C200, using the formalism of iMo. Mao fe White! 
(1998). Adopting an observationally motivated local mass- 
to- light ratio, (Af/L)tot, we identified those haloes most 
likely to contain LSB galaxy disks, thus allowing us to 
construct a distribution of concentrations that we compared 
with the distribution of inferred concentrations presented 
in DB03. Although there are important differences between 
the DB03 and our hypothetical LSB distributions, namely 
the presence of a significant number of galaxies for which 
C200 = and several high-C2oo systems, we recover a peak 
value of C200 ~ 6 that is in good agreement with observa- 
tions. In addition, when we exclude the C200 = systems, 
which are not true measurements of halo concentration 
but rather a failing of the entire fitting formula, and the 
C200 ~ 19 system, which likely has a much lower mass 
than our haloes, we find that the observed and predicted 
distributions of C200 are consistent at the 15% level. 

This is a compelling result because it establishes the 
relation between a halo's concentration and its spin, and 
confirms the expectation that haloes with low C200 and high 
A are the natural hosts of LSB galaxies. This is in good 
agreement with the observational evidence, which indicates 
that LSB galaxies occupy dark matter haloes that are more 
diffuse and more exte nded than those of HSB galaxies of 
comparable mass fe.g. lBothun et alJll997h . As noted in § 
this is equivalent to stating that the host haloes of LSB 
galaxies are less centrally concentrated than those of their 
HSB counterparts. B01 observed and quantified a scatter in 
concentration for fixed halo mass, and found a 1-er varia- 
tion of Alogcw ~ 0.14 (although we note that this has 
been revised to AlogCi,i r ~ 0.18; Bullock, private com- 
munication). Thus we expect a galaxy halo with a virial 
mass M200 = 2 x 1O 12 /i _1 M0 to have a concentration in 
the range 6 ^ C200 ^ 12.5 (B01, ENS). However, previous 
studies failed to identify any correlation between a halo's 
concentration and its spin; NFW examined the distribution 
of residuals in characteristic overdensity 8 C (a function of the 
concentration ) and fou nd no correlation with spin. Similarly, 
Bullo ck et al ] fcOOlbft found no evidence for a correlation 
between the spin parameter A and concentration, although 
they noted a weak trend between their modified spin A' (see 
equation|3J and concentration c v i r , such that smaller c v ir cor- 
responded to larger A' and vice versa. We suggest th at the 
absence of a correlation in the iBullock et alJ 1200 lbl) data 
reflects the relatively poor mass resolution of their simu- 
lation in comparison with ours (a factor of ~ 10 in mass), 
while we suspect that the selective nature 2 and poor number 
statistics of the NFW study will have made the correlation 
difficult to identify. 

Although this result does not resolve the current de- 
bate surrounding the mass profiles of LSB galaxy haloes 
and the implications for the CDM model, it is promising 
that a set of reasonable assumptions and a simple correla- 
tion between concentration and spin can so elegantly lead 



2 NFW studied a biased sample - systems in dynamical equilib- 
rium. 



to two halo populations that bear the characteristics of LSB 
and HSB galaxies. The central density profiles of the haloes 
must still be reconciled with the detailed shapes of LSB rota- 
tion curves. One possibility is that an additional correlation 
exists, such that LSB galaxies occupy I0W-C200, high-A and 
triaxial haloes, which may have observational consequences 
for the kind of rotation curves that we might observe; indeed, 
the bearing of triaxi ality on measured rota tion curves has 
been investigated bv lHavashi et alJ a2004ol . Intriguingly, a 
subpopulation of our LSB host haloes have particularly tri- 
axial figures. Overall, we consider this an encouraging result 
for advocates of the CDM model and one worthy of further 
study. 
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